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AgSTI-fjjgJ! 


Streariline curvature curd the matrix through flow 
methocic autempt to model s, real flo’iu in a turbomachine 
by a combina.tiOi: of .an axicyrmietr icj im-rscid flow 
model and a cascade lore model, ihe success of flow 
predictions depend on tfie adequacy of the loss model. 
The present ' orb: describes an application of stream- 
line curvature method to predict the perfonrxance of 
highly loaded multi- and single stage turbines using 
different loss model, A loss model has been proposed 
in this vYorlc which not only evaluates the overall loss 
level in a, turbine but 8,lso distributes these losses 
radially in the turbine annulus in a proper manner. 

The predicted results using this loss model 
compare very well itn tne experimental measurements 
for the above turbines. Results obtained using other 
loss models have also been comps.red \/ith those of 
present loss model. 
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CHAPTER 1 


IITROB UGTION 

Phe xlrst step in the prediction of the perforniance 
of a t art oma chine is to analyse the flow in it, (The internal 
flow in a turhoma chine is, in general, three dimensional, 
time dependent, and viscous and thus extremely complex to 
analyse. In addition, the flow is affected by such factors 
as tip clearaaice, coolant injection (in case of turhlnes) 
and by irregularities in the machine boundaries. 

Due to the conplex na,ture of the flov/, a solution of 
the flow equations taking most of the above flov? features 
into account simultaneously is a formidable task, Efforts 
at obtaining an analytical solution to the flow problem, in 
the past, have been made by making approximations that 
simplify the flow equations and permit the use of two^- 
dimensiona-l techniques for the solution. Some of these 
approximatiuns are i 

(1) The flow is assumed to be inviscid ai'id steady 

(2) Galciilation is confined to those regions of the 
turbomacnine where the axisyirmetry assumption is likely 
to be valid, such as 'duct flows' between blade rows. 

(3) Flow in the blade rows is taken to be circumferentially 
periodic with a period equal to the blade spacing. 
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(4) In analysing, the blade~t o-hlade flov/ in a rotating 
blade row, a frame of reference fixed to the rotating 
blade it cboaen and the flow is assumed steady with 
respect to it. 

(5) Boundary/' layers are dealt with separately and their 
effect on available flow area, flov/ deflection, and 
irreversible losses are iicorporated into an otherwise 
iiviscid calcula.tion. 

As first shov^m bj" (l), the t liree-dimensional 
inviscid flow in a turbomachine can be separated into two 
simpler but interrelated two-dimensional flows. One is the 
* blade-to-blade ’ flow v/hich occurs on a stream surface 
containing the tangential direction and describes the turning 
effect of the blade. The other is the meridional flow which 
occurs on a stream surface contaming locally the radial 
and axial flow directions and describing the radial equilib- 
rium of the flow. I'igure 1 gives a representation of the 
blade-to-blade and meridional flows. This splitting of the 
three-dimensional flow into two t'.io— dimensional flows has 
been fundamental to the success of turbomachnie performance 
calcula.t ions . 

The oresent work is concerned with the solution of the 
meridional flow problem to predict the performance ox an 
axial flow turbine. 
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1 • 1 Metli ods^ for so lyiiig ^ £2£>sL f ro t>lgj ii 
1.1.1 ^La-:^r ix Throug li glo v/ U© 


A_ 


Tae inethod involves tbe solution for the stream 
functioai on a curviJ-iJiear grid on the meridional plane, Ihe 
technique involves covering the region of interest v/ith a 
fixed irregular grid and vriting a finite difference approxi- 
mation to the principal equation, ^iven helov\r, at every 
interior grid point. 



g.(x,y, 


M. M 

6x * &y 


) 


( 1 . 1 ) 


This \/ill result in an algebraic equation for every 
interior grid po.int in terms of the stream function at that 
and neighbouring points. This system of equations can be 
expressed iii matrix form as 


Ca]] CO = Lsj (1-2) 

w here CaC 1-- the coefficient matri^c derived by replacing 

O - 

the differenticvl operator V ( ), L p J is the vector of 
unloiovTO stream function values, and QqJ is tne vector of 
quantities q(x,y) from equo.tion (l.l) and the boundary values. 

Since the right hand side of equation ( 1 . 2 ) is a 
function of 14^ and its derivatives, the system of equations 
is non-linear and must be solved iteratively, by first 
estimating LO , computing CO. and then repeatedly 
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solving etiuation (1.2) for ^ . The value of [Iq] 
improved in each iteration using the previous value of elf'll * 

This method of solution was first implemented by 
Marsh (2), Biniaris (3) has improved upon the formulation 
of the eq^uations by incorporating the concept of a 'represen- 
tative stream surface* and solving for the stream function 
on this surface. The blade forces v/ere determined as that 
required to Ireep the flov/ in the surface. This method 
improves the solution for blades with lean, 

Katsanis (4) has used the matrix method for obtaining 
subsonic solutions on a mid-channel surface in axial and 
m ixe d flow t urb oma ch ine s . 

Bosman and Marsh (5) have made further modifications 
to the formulation of the problem by v/riting the governing 
equation in a direction normal to the plane containing the 
velocity and body force vectors. By this, they are able to 
eliminate the bodj^r force terms from the equation. Furthermore, 
Bosman and liars h ( 5 ) have introduced a dissipative force term 
in the governing equation to remove an inconsistency in the 
formulation by the introduction of a loss model, 

1*1*2 The Fi n ite Jglonient Xie th od 

Tile formulation of the equations in this method 
is the same as that of the Matrix method. However, Instead o. 
using finite-differencing techniques as in the Matrix method, 
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the J’in.ite Eleuient methofi is used in solving the governing 
equations. Hirsch and harzeo (6) and liirsch (7) have 
described the method in detail. The riietaod is comparatively 
new and in the process oi development arid therefore yet to 
establish itself in turbomachine performance calculations. 

1.1.3 The S treamlixie C urvature _M gt lpod_ j S CTll, 

In the strea inline curvature method, the radial 
component of the If avier-S totes equations is expressed in terms 
of the fluid properties and streamline geometry and solved 
in a specified direction, either along a quasi-orthogonal 
or a radial line, for the meridional velocity. The form of 
the equation solved for an arbitrary quasi-orthogonal 
direction s, is given by : 



ds 



(1.3) 


whore 0^ and G2 a-^e coefficients 
orthogonal direction, s, and the 


which are functions of the 
meridional velocity, 


Since ecuation (l.3) requires a constant of integration, 
it must be solved in conjunction with the continuity equation, 
for the velocity along specified quas i-ortliogonals , 


Thus, there are bvo levels of iteration in the streamline 

curvature method, the outer iteration changes the streamline 
position based on the velocity distribution at each station, 
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and ttie inner iteration computea a velocity distribution 
which satisfies the continuity condition at each specified 
quasi— orthogonal. 

Smith (S) has derived the exact radial equilibrium 
equation in s, 3,’orni suitable for tne streamline curvature method. 
The equation has been derived by circumferential averaging 
of the flow variables in the moment um equations. 

Novaic (9) has described in detail the process of solving 
the meridional flo?/ problem using the streamline curvature 
method. He has derived the governing radial equilrbrium 
equation from the radial momentum equation of an inviscid 
flow by making use of the axisyiionetry approximation. Novak 
has also discussed the problem of a non-axisyiimietric 
computation which is necessary when computat ions are to be 
made within a blade row where the axisymsetry assumption is 
not valid. V/hen such a computation is attempted, terms 
involvkig derivatives in the circumferent ial direction 
occur in the equations and the handling of these terms poses 
f ormidable problems . 

The streamline curvature method for the off-design 
analysis of axial flow compressors has been implemented by 
Jansen and Koffat (l0)» The estimation of losses was 
accomplished by using cascade correlations developed by 
Mellor and Liebliti (li). Corrections have also been made for 
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the effects of streaBilijie slope, axial velocity variations, 
"bla,!© tbiclcnesB, and Macli numbers on the correlations - The 
paper is a pointer to the use and suitability of the stream- 
line curvature method in turbomachine computations. 

Frost (i2) has improved upon tue metbod so as to be able 
to make computations within blade rows. This has been 
achieved by using the concept of a 'raean stream surface’ 
developed for bhe Matrix method. The flow equations have 
been derived on this stream surface and the blade forces 
determined as tliat required to keep the flo'w in the surface. 
Frost has called this method the ’streamline curvature through- 
flow method’ , Losses have been accounted for by prescribiiig 
effioioncles for each blade row. 

In a more sophisticated analysis, hovak and Bearsley (l3) 
have obtained a ’nearly three-diiriensional’ solution to the 
flov/ in a turboraachine by using the streamline curvature method 
and iterating bet\'/een the solutions obtained on the blade- 
to-blade aaid meridional pianos. Computations can be carried 
out within blade rov/s ’^’itii lean .md twist, Hovi/ever tjie 
analysis is restricted tc the invisoid flow. 

An application of sai to multistage steam turbines was 
made by Renaudin and Soram (i4). The method uses the blade 
height for a given mans flow as the boundary condition 
instead of the mass flow for a given blade height more 
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commonly used. Slneamlines have been assumed to he elastic 
liases 'vhich are a function of time instead of space used in 
similar spline fitting techniqueB. 

Katsaiiis (4) has used the method in conjunction with 
the i^atrix method for evaluation of transonic flows ivhich is 
not possible using the Mafcrix method. Silvester and 
Hetherington (l5) have also developed the method independently 
on the same lines as that of liovalc (9). 

In an application to gas turbine design, Carter, 

Platt and Lenherr (l6) have used the SOI for the evaluation 
of the design point perf ormancc of axial turbines and also 
as a tool to evaluate different des igns .during the design 
stage of the turbines, ilie radial equilibrium equation has 
been simplified by neglecting the term involving the change 
in momentum in the meridional direction by assuming that the 
stixsamline slopes are small, ihe slope and curvature of 
the streamlines Y>fere evaluated from the slopes and curvatures 
of the hub and ca.sin:; boundaries: by assuming a linear 
variation frO”: the hub to casing for them. Ihis is in 
contrast to the ” spline- fit ting” tedmique adopted by 
bovak (9), Jansen and 11 off at (lo), Silvester and he the ring ton 
(15), among others. A loss correlation has also been 
developed by evaluating the pressure loss coefficients from 
the blade angles of the blade row and correlating thori with 
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tile efficienoy plots of Smttli (IV). Effect of iiiterf ilaiiient 
sii^iTig aii'i cooltant injection lias also been included in. bbe 
Gomp ater prograa , 

1 . 1 . d gy:^is orj_ p f _b h®. aml ine Obryaturg _ Meth ods 

A co'uparison of the matrix method ana the streamline 
cur-vature method has been made by Davis and Millar (18). One 
of the important conclusions dravm by them is that it is not 
possiOle for oitner motbod to claim a clear operational 
advantage over the other. Some of the advantages claimed 
for the matrix method arc slightly better accuracy aiid 
stability of the computing procedure, Hov^ever the matrix 
method requires a large computer storage for implementation. 

Some of the advaiitages of the streamline curvature method 
are its ability to handle transonic flov;s and small computer 
storage for imp 1 orientation. Another important conclusion 
drawn by Davis and Millar has been that tho accuracy with 
>7biGh tho xlo\'7 is predicted is ultimately dependent on the 
cascade model used for tho evaluation of losses in a blade 
r ov; . 

liie ability to handle transonic flovjs is an 
important reduiroment for axial turbine analysis as the 
flov/ at tho first stator exit is invariably transonic. Along 
v;ith this, the need for only small computer storage require- 
ment led to the choice of the streamline curvature method 
in tho present v/ori. 
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1 . 2 Los 0 m odel 

The eq.uationB ox I'lovr used in the meiliods discussed, 
above ^ are for an invlscid flor/, Ho\';ever, t>ie analysis of the 
inYiscid flovr would be of little use in. the pradictioii of 
performance of the turboiuach ine unless some method is found 
to take into ancoLmt the clianyes in entrop/ due to Yiscosity. 
This is the ouroose of the loss model. The loss model 
superooses on the inviscid flow the effects oi viscosity 
in the form of losses occur ing in the flow. These losses 
may he superposed in the form of a stagiiation pressure loss, 
an increase in entropy, loss in kinetic energy, or hy the use 
of efficiencies. 

The evaluat ion of the losses in the flow is usually 
achieved from correlations of experimental data. A large 
amount of such data is amtiable in the case of axial 
compressors. A NASA special publication (i)ASA SP~36) on 
axial compresooro contains (hi. collection of such data. 

Masumdar (l9) uas revxe^ved the different correlations 
available for ourial compressors. In the case of axial 
turbines, however, loss correlations are limited. Tv/o of 
the wellknoYm and much used correlations are the Ainley- 
llathieson correlation and the Soderberg correlation 
(Horlock (20)). 

Another method of evaluating the losses in the flow 
is by solving the boundary layer eqtUations. 


This is the 
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current trend in turbomachine calculations, Tbe end-v/all 
boundary layerc are three-dinene xonal ana the solution of the 
fare e-diiaeneioiial boundary layv,-r equations is too complex 
for impleraentat ion into a design- analysis computing system. 

While the evaluation of losses is in itself importantj 
the distribution of these losses along the blade height has 
also to be tfnen into account. The distribution of these 
losses must be determmed from experimental data. Unfortunately, 
such data is not available in published literature. The 
radial distribution of these losses, therefore, should be 
made in a simplified manner consistent nith the nature of 
the florv. Jansen (2l) has made each an attempt for cejitri- 
fugal compressors. The skin friction losses have been 
distribiitea in proportion to the \7etted area and square of 
the velocity. The clearance losses have been distributed 
befveeii the first four streamtubc->s , at the blade tip, in a 
specified proportion, hesulto obtained hov/ever have not been 
compared against experimental .lata. 

P i G sent inyp s tiga tiong. 

The streamluifc: curvature method is a well established 
method for turbomachine flov' analysis "»/hich can be implemented 
on a rolat Lvely small computer system. This along "v/itb, the 
fact that the method can be used for transonic flows led 
to its choice in this worl: for the design-Doint performance 
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prediction or axial turPxnes. 

She accuracy v/ith which the performance can be 
predicted deoends to a large extent on the loss model used. 
As already explained, it is also necessary that these losses 
be distributed appropriately along the blade height. A loss 
model using 7'ell icnoiTO loss correlations but also taking 
into account the distribution of these losses has sou,,ht to 
be deyeloped here. 


I 




jjg^.IIJLATIOB ^PP^‘ACxI_ TO^ SOIil^ M 
2 , 1 gormul-_t lon 

The inviscicl tliree'-diraensioiial flovr pattexn can be 
found by sol'/tn/'^, the follo^/ung equations 

( l) O-U'itinulty (i equation) 

(S) I'.'iotioii (3 equations) 

(3) Energy (l equation) 

(4) State (l equation) 

Solution of these siz equations determine the fluid 
properties f , h and s and the three velocity componentG 
V , V and T . For a reversible adiabatic flav the three 

-C W ,A. 

equatio 2 is of motion can be combined v/ith the energy equation 
to obtain an equation that states that the entropy of the 
fluid remains constant along any streamline. Since the 
equations of notion, energy and entropy/" are interdependent, 
solutions can be obtained by using aiiy four of them, the 
fifth being automatically satisfied. 

If the reversible aliabatic flow takes place on a 
prescribed stream surface (the meridional surfai^e), then 
there ezlste a geometrical condition relating the three 
components of the velocity vector in order that the velocity 
vector lies in the stream surface. 



Vfitli this sixth equation, there is a roduiidanc.y ot 
one equation and i^ne flow pattern is oTofcained hy replacing 
one of the equations of motron with the geometrical condition 
for the to remain in. tae stream surface- fhe replaced 

equation of motion can he used to evaluate the force necessary 
to keep the flovif in. the stream surface. This force, for an 
inviscid flow, lies normal to the surface and is contained 
only in the equation of motion tor this direction. Thus 
this equation is replaced hy the geometrical condition for 
the flo^./ to remain in the stream surface and can he used to 
evaluate the force required to do so. 

Thus for a reversible 9 .dtahatic flow on a preset i-ood 
stream surface the flow pattern can be determined by the 
following six equations 

(t) Cont;inuity 

(2) (a) iiny two of ; Zi equation of luotion 

r e 'V uat i. on of m o t io ii 
mntropy 

( b ) G-e omet r ioul coiiv. it i on 

(3) Energy 

(4) State 

S’ or an irreversible adiabatic flow, tho entropy 
increases along the streamline and this effect is talcen into 
account by means of the lose model* Since the irrevers ibilrt y 
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of the flow is caused hy viscous effects^ the ijitroduct ion 
of a loss model superposes the viscous effects of the real 
flow on m\ otherwise Inviscid -formu] ,at ion, 'This leads to 
a certa.ln iai consistency in the equations (Horlock; (Tii))as 
the equa.tions are "basically for an invir. cid flow. An 
alternate foritulatijn to avoid this inconsistency, proposed 
hy Hoxloch (2t) and also used by Bosnian and Marsh (5) has 
been shown to he invalid hy Wu (23), Also, in spite of this 
inconsistency, tne numerical solutions for the flow pattern 
by using loss models have been very successful in 
meridiona.1 flow predictions. It is this set of equations 
teat are widely used in the anolysis of the flow in 
turbomaxhines and has been used in the present investigation' 

Thus, for an irreveisiblc adiabo.tic flow on a 
prescribed stream surface, the six governing equations are j 

1) Continuity 

2) a) l.iotion (r direction) 

b) Entrooy (loss model) 

c) G-eometrical condition 

3 ) Ene rgy 

4) State. 

The equation of motion in the r direction describes 
the equilibrium of the floir along tne radius and is the 
y/ell Imown radial equilibrium equation. 
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2 . S ion s 

2 .<■ . 1 iuiTi^ ej^uat ion 

T^'- ILe CiOrivat. ion oi tlso radial- oquilil'iiuii e^^^ual ion 
the iollowiiife as sunipt ions have been Jiiade . 


1 . 

2, 


r5 


4. 


The now in. hiviscid. 


T he rot or 


is rip. id and rotates with constant 


angular 


velocity 


T]je llnid is a seiaiper.fect gas x.e. the equation of 
state is P =: pHT (2.l) 

vjitb H Gonstaait, and the speci.rio beats are 
deoendent upon the temperature only. 

The flov-; is steady and a:jnsymmetr i o , 


Por a frame of reference rotating v/ith constant 
angular velocity ui about the aris, Eowton's Becoiid lev/ 
of motion gives for sii inviscld fluid, (Pigure (g)) 


VP _ in 

p " Dt 


2-^ 

w r + 


2w X 


If 


( 2 . 2 ) 


The unit vector derivatr-res are 




and 




i 

r 


and the hhieniatic relation 2^^- = . 

Ut Ds 

The .0^' -operator signifies that a fluid particle is 
being followed during different iat ion. The equation (2.2) 
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call be writ cell iii terms of its scalar comiorentB as. 


1 ^ 

'p 5? ' 

1 

p rde 


('A'n + wr)' 


•Oi/ W 

Tj ^ u r o -T 
~ + 2a, 


i ^ „ IV ^ 
p 5i ' Ds“' 


(2,2a) 


(2.3b) 


(2.5c) 


l;3ir.£i tile substitution :=r -f tor betiveen tjie absolute and 
relatiTe velocities, equations (2.3a) and (2. 3b) can be 
rewritten as. 

1 6p , 

- V 15 = Vs~ - ~ 

_ 1 1 i£ - D , 

p r &e J \) (2,4b) 


We can als o ~i r it e , 


I?/ ?XJL _ ly 
■ Ds " X Dx 


(2.5) 


where Dx represents the increago in the x coordinate tVie 
particle undergoes as it noves o, lis canoe Ds in the flow 
direction. iJiua-tions like (2.5) can be v/ritten for aiiy 
desired, direction. Written for the meridional direction, 
defined b.y 


i Dm = i bx + Dr 
m X r 


(2.6) 


equat ion (2.5) he comes 

1.. I _ i/ ?i-i 

■' Bs " 'm Bm 


(E.7) 
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Subat itucing equation (2,7) in equation (2.4a), we get, 


. . Y 

1 op ^ __r _ _U 

p or ' TJ- Dpi r 


further, we can write =: oin 

and equo,tiori (2,3) can be rewritten aj 


( 2 , 8 ) 


(2,9) 


1 PZ ^ . V- 

p or r , 


D slutj) 

Dui n Dill 


( 2 . 10 ) 


D Gin ^ 


- =: COG 4> 


( 2 , 11 ) 


'Due radius of curvature is considered negative when the 
streamline proaection on the meridional plane is concave 
up^jond , 

aquation (2,io) can no/; be rewritten as 


^ D7,. 

1 ^ 'C) u , m -j ^ 

i et = r' cos* - q 555- 


(2.12) 


V/ritten in. terms of the ahsoluiio velocity coHiponents, 
e q uat i on ( 2 .12) be cocic s 

Y^ D?.., 

1 i£ ^ ^ 3 „ V sin 4, (2.13) 

p dr r 

Bq uat ion (2.l3) is written in a more convenient form by 
replacing the stottc pressure by the stagnation pressure from 

„ yS y2 y/y-1 
p = Po ^ " TS*^ 

!— p O ■"* 
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Je get 


U iii 


1 <3.2 dv,. 

■1^. -P- ^ Y - 

5q dr " u dr 


ittS 

-1 

1 m 

C' 

KJ 11 X 


+ E T, 


o 


p p 

Y +Y 
u m 

Sc 1 . 

p o 


X 

P. 


dp 

dr 


0 


2 ? 

!iu dli 

r r 

1-1 


y. 


sm(j) 


Dy. 

1 

Hr 


ili 


(2.14) 


m 


It 13 HOW necessary' to express tne term in terms 


of the raeridioiial pathline slope and curvature, dy making 
use of the continuity ec[uationj erj_uation of state j and laws 
of tbermodynamicSj it can be shown that (Appendix A) 


DY.. + % + j. 




sin 4 > 


m 


jyooF?^ ♦ aJD 


Em 


1 - k, 


hi 


(2.15) 


Smith (8) has derived the above e liquations by clefij'j.ing 
averages of the flov/ properties in the circumferential diiectioi. , 
It has also been shovm that the enuations are the same as that 
de r ive d b y h ov ah ( 9 ) b / mol: j i\: x he ax is yuuie t r ic ap p r ox imat ion 
when the ” G-~ functions” ■1e-finecl by Smith (8) can be neglected. 

It in evident iror., tue analysis of Smith (8) that the evalua- 
tion of the defunct ions is too difficult to be incorporated in 
a meridional flow solution method. Another difference between 
the approaches of Smith (8) and hovalc (9) is in the choice of 
the dependent variable in the radial equilibrium equation. 

V/hile Sraith has used the static pisssure, 11 ovale has used the 
meridional velocity as the primary dependent variable. This 
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Is purely a matter of convenience and the meridional velocity; 
as s.dopted hy lTovs.k (9), has been used here as the dependent 
var iahle . 

2.2,2 T he Go ntiiiLi lty Bqua txon 

The continuity equation provides the necessary condition 
for the deteriflination of the constant of integration when 
integrating the radial eq_u 11 ibr ium equation. In other words 
it provides the necessary boundary condition for the radial 
equilibrium eoy’.ation. Thus while the radial equilibrium 
equation determines the distribution of the meridional velocity 
across the annulus, the continuity equation determines the 
magnitude of these velocities so as to satisfy design mass 
flow specifications. 

for an axisymmetric flow passing through an axial 
section of an arbitrary aimulus, the continuity equation can 
be written in its integral form as 

^t 

W = 271 j" p T. r dr (2.l6) 

IL ji. 

The equation is v^riLton in a more convenient form in 
terms of the to'hal temperature, total pressure and tne 
tangential and meridional components of the absolute velocity 
which hove been chosen as the principal variables in the 
analysis. The density can be ’written in terms of the total 
pressure and total temperature as 
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o 



T® -1 
pc-*'* 


(S.17) 


SuDstitut© (2ti7) in (S»i6) to 


_ 2 tz pQ 

Wj - J — 
-h ^ 


+ V^’ 

U lil 

2 c t 
P 0 


Y-1 


V,., cos <j> r dr 

ui 


(2.18) 


2.2,3 Tbe los s egu at io^i 

Since t’le total teraperf ture and total pressure are 
chosen as toe principal variables in the analysis, the 
irreversibilitT of the flov/ is more conveniently expressed 
in terms of a, sta^ynation pressure loss than in terms of the 
increase in entropy, For this purpose a total-pres cure-loss 
coefficient is defined as follov/s ; 


or 





1 + 




(2.19) 

( 2 . 20 ) 


For a rotor, the loss coefficient is deiu^cd as (FioUX'e 3) 


- 


o2; 


■R 


n 




- p 

-“pT; 


o2 


( 2 . 21 ) 


or 



1 + Y^(l 


p. 


o2 


( 2 . 2 ?) 


Ihe above eruaiions can be used to determine the stagnation 
pressure loss \7hen the totaL-presaure-^loss coefficient is known. 
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The loss coefficient is CTOvluated. from the correlations of 
experimental loss data in the form of a loss model. 

If toe optioi- of see c if p- me; stare isentrooic or rotor 
isentropic eff ici-enoieG is used, equations similar to (2,2o) 
and (2,28) can he v’/ritten for toe etgige outlet stagnation 
pressure 


■^ o 2 -^00 


1 


f f ■ 

“77 

f S 00 


p 


o2 ^ol 




- f 

Q l Pg 
1 f 
t R 00 


1 

y~T 


Y-l 


(2.23) 


(2.24) 


2.2.5 Ttie \/p rli e quation 

The necessary equation rela/bing the tangential 
Telocity to other principal yariahles is provided, at a rotor 
er.it, by the Suler vork equation. 


or 


A = V,,. u, - u. 


ul 1 


2 


V 


u£ 


- c AT 

Ul p o 


(2.25) 

(2.26) 


This eouation is used in the analysis v/hen the design 


oil: is spe cif is d . 


2.2.6 pond i t i p ips 

At a stator exit the tangential velocity is related 
to the liieriili jnal velocity from a specification of the exit 
flow angle in the Telocity triangle. 



23 


Y = Y cos ^ tail 3 (2,27) 

U. LI 

Vhen the coTaplete '‘relocity triaii-tles are specified, 
the taiigent ial velocity at a rotor exit can be related to the 
meridional velocity by oji equation similar to (2,27). This 
equation replaces the Surer ivorS equation for the tangential 
velocity in tae aneSysis ot the flon » The Suler v/ork equation 
(2.2 5) can no\; be used ior the evaluation of the total 
temperature Iroq through a stage and hence the total temperature 
at tue stage exit . 

2 , 3 Appr oach to the solution of t he e quations . 

The solution of the flow field in a turbine conpisto 
essentially of obtaining values of the total pressure, total 
t e up c ref u re } and the t',/o components of the absolute velocixy 
at each selected streamline location in each of the radial 
planes e.t vmich the solution is being obtained (design plane). 
Thus, Cor tie onalysis p^, and Y^, are the principal 

variables. In the previous section all the equations Yi.iich 
have to be satisfied in the analysis have been presented In 
terms of these variables. The teciinique for solution adopted 
here is similar to that by Carter, et , al. (l6). 

The set of equaiions, for the solution, would comprise 
of equation (2,li), tae appropriate equation in the set of 
eqr.ations (2.19) to (2.24), equation (2.25) or (2.26) and 
equation (2.27). It is obvious that the principal equation in 
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t h is set is t bs rad ial ec^u il ib i ium e quat ion ( 2 . i4 ) , T li is 
equation is aai ordinary differential equation, nonlinoar in 
the varia,tlo The remaining of the equations can xio\7 he 

used for the evaluation of the other flow variahleo in equation 
(2.l4), The total temperature is evaluated from the Euler vrork 
equation (2,25). (Appendiv B). In addition, the streamline 
slcoe and curvature in the meridional plane ne cessita-tes 
consideration of derivatives with respect to the axial direction 
X. That is s 


^ = tan“l 0 




,2 2 


d^r/ di 


r. 




Hid- (dr/dx)^']^'^ 


(2.28) 

(2,29) 


Kcnce, the radial equilibrium equation (2.i4) contains 
derivatives with respect to hots r oiid x. It is unrealistic to 
assume that the axisymuetr ic form of the radial equilihrium 
equation can he extended heyond the inteihlade ro'W space liito 
the hlade ro'vs. Thus, the boundary conditions for the 
meridional streaitilin 3 b hi the interhlade row space are indeter- 
minate at the trailing edge end leading edge planes, defining 
this apace. Only the houndary streamlines, at the inner 
and outer W3.11s, are defined hy thi assumption that these 
streamlines follow the contours of the annulus wall* In the 
absence of a rigorous analytical treatment for the slope and 
cur’/ature of the streamlines in the meridional plane, it 
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becomes neeessai’y to make more assumptions for the solution of 
the problem. In the present analysis, the streamline slope 
and curvature in tJie meridional plane has been calculated in 
one of two ways. 

a) The slope end curvature of the boundary streojiilines 

be obtanied from a definition of the wall contours, and 
then Doth ^ and r,.^ will be assumed to he linear 
fuiictions of radius determined from the values at the 
v/ai-ls . 

b) The use of spline- fit ting techniques. 

Detailed discussion on these methods is presented in a later 
chapter. 

A nunorical solution of the equation (2,l4), clearly, 
has to be iterative. The continuity equation (2,i6) provides 
the necessary boundary condition for the evaluation of the 
constant of inteyration in equation (2.11). 

2.4 The djL^:^r^8 irfc ia 1 ejqu^t ions 

In precentxiyj tue equations here, the case of a single 
stage turbii.e has b.jen considered v/ith design stations at the 
stator inlet, stator exit, and the stage exit, A faultistage 
turbine can he considered as mereljr repeating tne solution 
technique of the first stator exit aiid the first sta-ge exit 
planes . 

The radial equilibrium equation (2,l4) can be written 
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in tile fora 


n 4, n i--- -=-=-^ 4. n — - 

11 dr ^ ^12 dr ^ ^13 dr 


(2.30) 


Y/liere the coefficient C^g, 0^^ end 0^^ can be assigned 

values at each point in a particular design plane once a value 
of neridional velocity is selected. Ihe variables T^, r^ 
and r and constants Y and are assumed to have been evaluated, 

The differentiation of the appropriate equation for 
total pressure (Any of the equations (2.2o), (2.22), (2.23) or 
(2.24), depending on the station and option being considered) 


leads to a differential equation also of the form 

p 

^ "iTi* 


'21 dr 


+ C, 


1 


22 p dr 
^ o 


+ 0. 


23 dr 


(2,31) 


where the coefficients ^22’ ^23 ^24 evaluated 

in term'? of the input specifications and meridional velocity. 

Similarlyj differentiat ion of either of the equations 
(2.26) or (2,27) for results in a differential equation 
again of the form 


n 4. n -1^ 4. n - n 

^31 dr **■ ^32 p dr ^33 dr " ^34 

o 

In this case coefficient Cgg is alv/ays zero. 


(2.32) 


Ihe terras involved in tlie coefficients are presented 
in Apoeiidix (C) after differentiation of appropriate equations. 
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2 * 5 Techn ique f or solut ion 

Ihe re’vrit mg of toe oqaations of Chapter 2 in the 
form (2.30), (2,3l) and. (2 ,32) helps ii: develop ing aii unique 
technique for the solution irrespective o.f the design station 
beini considered, only the coefficients ol the dlfforont ial 
equation differ from station to sta.tion. the problem of 
solution is Jiov/ one of ohtutnina o. meridional velocity 
distribution which simultajioously satisfies the radial 
equilibrium and continuity enuations. By choosing on initial 
value of the meridional velocity at one strearalinc position, 
the local values of the coefficients of the set of equations 
( 2 . 30)5 (2.31) and (2.32) can be obtained. These equations 
are nm solved for the derivative dV^/dr. bsing the standard 

ui 

Runge-I'Cutia-Gill method of solviiig ordinary different in 3. 
equations (Appendix D) , the value of the meridional velocity 
at an adja.cent strecuml iixG can bo obtained. (Jsing the tie\: value 
of the Tiler id ional velocity, values of the to Lai pressure, 
tangential velocity’-, and tue coefficients can be obtained 
at the new streo-mline . Tjus the dc.rivatlve d? /dr can. be 
obtained at this nev streamline also. The process is repeated 
until the iSo L’idional velocity, total pressure, and tangential 
velocity are obtained at each of the streamliiios used in the 
analysis. Using the continuity equation (2. 18 ), the mass 
flow for the meridional velocity distribution obtained is 
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evaluated. Since the oolution was obtained froii an assumed 
va.lue of the meridional velocity at one point in the flow 
field, the mass flow computed ill in general differ from 
the design specified value. Hence, the assuiaed value of the 
meridional velocity 'will have to be modified iteratively until 
the starting value is consistent \rith the continuity require- 
ment, This constitutes the inner loop of the iteration 
procedure at a design station. 

The outer loop of the iteration scheme is dependent 
on the method adopted for evaluating the streamline slope and 
curvature. First, we assume that method (a) of section 2.3) 
above, is chosen for this purpose. In this case, it is clear 
that the slopes and curva.tures are functions of only radius 
at a given design station. It is now possible to complete 
the outer loop of the iteration scheme at the design station 
before proceeding to the next station downstreani, This is 
completed as follows i- Since initially the flow distribution 
is unlmovTO, the initial streamline positions are estimated 
from eoual areas for each stream tube. Ilenco, streamline 
positions have to be relocated after each solution of the 
radial equilibrium and continuity equations until a converged 
solution for streamline positions has been obtained. This 
completes the solution at a design station and the analysis 
can proceed to the next design station dovmstrearn. 
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In tiie case of adopting tbe iiiethod of spline fittings 
method (h) of section 8, 3, a nev/ iteration scheme has to be 
applied for tie outer loo), fh' roasoii for this is that in 
spline fitting, the slopes an;! curve^tures ere functions of the 
radial position of sxreanlines at upstream and downstream 
stations also. Thus it is naoegsary to complote the inner 
iteration loop at all design stations in the turbine annulus? 
the streamlines bo relocofced at these stations as in the 
previous methods, and a new set of spline ciurves be fitted in 
the annulus* Thus, after completing the inner loop at a 
design station, the analysis proceeds to the next design 
station downstream till the entire turbine annulus is traver- 
sed. Strea-mlines arc relocated at all the design stations and 
spline curves are fitted to simulato a new set of streamlines. 
The inner loon is again started at the first station and 
proceeds downstream till the annulus is traversed* This 
process is repeated till convergence of streamline positions 
is obtained. 



G HAPgGR 5 
LOSS 

3 . 1 ^ tj:pduct ion 

■riie tasL of predict uit? the design-point performance of 
of an axial turhine c£ir: he acco.upliehed only \t'aeii the viscous 
effects of the flov/ can be incoip or- ted into the basically 
inviscid calculation of the stroaiil me -curvature raethod* IhiSj 
as has been explained in Chapter is o.chieved by replacing 
one of the riiomentum equations by an equation which describes 
the loss of stagnation pressure in a blade ran. In using this 
equation (eqns. 2,19 or 2,23), it is clear that the pressure- 
loss-coefficient, Y, or the blade row efficiency, n , should 
be known beforehand. This, then is tho purpose of the loss 
model . 

It is self-evident at this stage that evon though the 
preparation of numerical procedures for the streamline- 
curvature method of performance prediction is of importance 
in itself, unless the procedure incorporates soundly based 
assumptions concerning the losses ascoGia,tpd v/ith the elements 
of the blading, the detailed calculan ions of interblade rovf 
aerodynamics v/ill be of little velue. Hence, an essential 
part of the development of a design analysis computiiig system 
is the development of a loss correlation which muld be an 
integral part of the computer program. 
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Unlike the case of axial compressors, loss data is 
limited in the case of axial turUiiies. Two of the most 
coLmonly used correlat iont are the Soderherg correlation and 
the Ainley— I'iath ieson correlftt ioii» These correlations have 
"been developed us ini; da,ta oUtained from steam turbines and 
early gas turhines. Uov/ever, these correla,t ions a,re still 
in use in modified forms because of lack of other reliable 
data (liorlock (20)). 

Three loss models have been used in this vnrk. T'vo 
of these have been developed by Garter, et , al, (l6) aa^id 
the third has been developed here to rectify some important 
drawbaclvS of the other t\70 models. Thevse loss models have 
been described in detail in next section of the chapter. 

3.3-1 lo ss model - 1 

This is the simplest of the loss models wherein loss 
data in the form of pressure-loss-coef ficient s , kinetrc- 
energy loss coefficients or blade offieloncies ore specified 
in the incut data to the comouter program, "dvidently, such 
a model can be used only \;heii data is available from experi- 
mental investigations on the turbine being analysed or from 
experience obtained from other similar turbine investigations. 
The use of such a loss model is restrictive and may be used 
in preliminary design analyses. 
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3 »3 .2 XiOoS model 2 


I’le loss correlation used to Uiis modol lias been 
developed by Garter j et . e.l. ( l6) . Ibe dotailed development 
of the corrolation is available m tnat reference. Tbe 
sta,rtin'5; point of tbis development is the Soderberg correlation 
gii/en belov/. 

i - (0.00076 + 0.00000 56^) ^ ^-9 ^ ^ 

* • » ( 3 • l) 

The above correlation implies th^xt v/lieii the velocity ratio, 

Y. /V , falls below 0.45, the level of loss for a given 
deflection begins to increase. This is contradictory to the 
expectation tnot the losses would decrease smoothly as the 
over-all row acceleration increases at a constant deflection. 

To overcoxic this discrepancy, Carter, et . al . Cl6) have 
assumed that the losses are proportional to the blade inlet 
and exit angles instead of oii the do fleet ion as in L-he 
Soderberg correlation. They pave developed a correlation of 
the form. 


tan p - tan ^ 


Y _ ra^+ao(Y /Y a„) > 




and Y = 


tan p - tan p 


(a. + Up. cos 


ex 


{ag + a^(Y^yY^^) ^ } 


if Y - /Y < 

ex 3 


(3.2b) 
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The 3asi& of the correlation are the efficiency plots of 

Smith (l?). Garter, et , al, (lb) have assumed that the 

efficiencies ohtainet5 from this plot are the achievable 

efficiencies of the desiyn being analysed. With this 

assumption, the constants a., to in the correlation have 

1 y 

been 6va.luated so that it prod tots the efi ic leixcies of 
Smitns (l7) plot. 

3.3.3 Lpjs 3^ mo del 5 


The loss model described here 1ms been developed to 
Overcome some importaiit shortcomings of the t\/o other loss 
models just described. As explained oarlier, the use of 
the first loss model is limited by the availability of data 
to use blade row efficiencies as inputs to the calculation 
scherae, X scrutiny of the development of the second loss 
model shons that the loss levels calculated from, the 
correlation are primarily dependent on the hlamle angles of 
the blade rov/ . Thus, the distribution of losses across the 
annulus would reflect the trend in the distribution of the 


blade angles. As aii example, if the blade angles are largest 
at tile hub and least at the tip, as is true in mgaiy turbines, 
the loss levels calculated us 11 % correlation (3,G) would shov/ 


minimum losses at the tip and maximum Icssos at the hub. 


This distribution is contrary to experimental observo.tions 


and qualitative studies based on the nature of lo-ses in a 
blade rov^ which show that the losses are minimum ai'ound the 
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centre or tlie annulus and incneaGO tov/a^'"^ ‘tti® tip 

(24, 25, 26). l-buG, vliila the correlation (3.2) ciay predict 
the overall oIlicioncioG and average Iogg levels accurately, 
it cannot predict the listr ihutioic of these locsos across the 
annulus . 

To analyse the a istr that ion of losses in the annulus 
Of a turoine, it is essential to seperato the losses into 
its prinoiual constituents, na-inely, the profile loss 'Und the 
secondary loss. The profile loss caused by the houndenj' 
layer on the suction and pressure surfaces of the blade is 
dependent on the blade an^^les (Horlock (2o) , Balje (25)). 

Thus, it v/ould be reasonable to assume that the profile loss 
calculated from correlations based on blade angles would 
reflect the correct trend of profile loss distribution in 
the turb ii'ic annulus. The secoiida.ry losses a,re caused by 
the t’nree-dunensional boundary layer at the hub and tip of 
the turbine annulus. Wh^ile it may be possible to calculate the 
overall level of secondary loss from correlations based on 
blade vt- their distribution in the turbine annulus 

cannot be obtained from such a correlation. Qualitative 
explanations and measurements (lianley (24), Balie (25), 

Hosney (26)) shov/ tnax the secondary looses are iiiazimum in 
the hub and tip regions and negligible in the centre of tho 


annulus . 
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Therefore in toe developnent of s. loss model , it is 
essential to calculate the profile and secondary losses 
separately and to tol:o into accoiint the difference in their 
distributions m the turoine annulus. In the prescait loss 
model, the profile loss ha.s iDecri ca^lculatod usin^ the v/ell 
hiioiMi Ainley-Iiathieson correlation (Horloch (2o)) 


Y = 
P 


a. 


p(P.=o) 


+ (Y 


df 







2 


(3.3) 


Ihe losses are evaluated at each streamline in the calculation 
scheme and since the correlation is primarily dependent on 
the blade ang,les, the distribution of the profile loss would 
be acGUiMto, 


Many correlations are a'ia.il; hie for the evaluation 
of S8cond0,ry losses (Dunjicum and Cajdo (27)). In their reviev/, 
Dunham and Camo (2?) hove concluded tha.t the best Gorrela,tion 
is of the form. 

Yg = b^b3^(cos ^g/co3 P j_) ^Gj^/b^/b^) cos^ pg/cos'^ 


function (6^/b ) (3.4) 

1 o 

She evaluation of the inlet boundary layer thiclrncss, 6^, is 
difficult in the present calculation scheme and the function 
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is assumed to be constant, Dunham and Game (28) have used 
a value of 0.0334 for this constant. 


= 0.0334 (■ 




cos ^ 


) 2 


(3.5) 


"n — '"1 

where Z is the Ainley-Mathieson loading parameter 

„ , \a 003 ' ^2 

^ = VTT'^ 'T'rT~ 

c cos 


(3.6) 


The secondary losses were distributed across the 
annulus in a parabolic form. 


Y (r) = A r^ + Br + C (3*'?) 

s 

The evaluation of A, B, and C requires three conditions 
-which haxe been assumed talcing into consideration the 
nature of distribution of the secondary loss in the turbine 


annulus . 


( 1 ) 

( 2 ) 

(3) 


The mass-averaged secondary loss obtained from (3,?) 
is that evaluated using the correlation (3.5) 

The secondary losses are gero at the mean streamline. 
The slope of the ll istr ibut ion is zero at the mean 


streamline . 


The three conditions determine the values of Aj B 
and 0 in (3.7) and the expression can be used to evaluate 
the secondary loss at any streamline in the annulus. Thus, 

the total loss at any streamline is, 

I.Yp + Yyr). 
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QHAPTEE 4 
QOiv-FU-rSR PEpGEAIvi 

~ 1 Qene ral D escript ion of t he^ C ompu ter ProgreJft 

Tbe basic equations v/Iiicli govern the design point 
performance of an axio.1 flon turbine and the uiethod of their 
solution has been presented in previous chapters. Ihe 
numerical method of solution is best implemented on a digital 
computer. The computer program for this purpose, used here, 
has been based on the computor program of Carter et , al, ( 29 ), 
A general description of the program is given belov;. 

The program has been written in Fortran IV for use 
on an IB!;I 7o44/i4oi data processing system. It is capable 
of analysing both single and multispool units (a maximum of 
three spools is allowed) and each spool may have upto eight 
stages. The absolute and relative ilov.’’ fields are computed 
at the first stator inlet, at each hitcrblade row plane, and 
at the final rotor exiit. The ehiects of blade coolant flows, 
interfilanent mixing and a stat ion-to-stat ion vaiiation of 
the specific heat can be included. The loss correlations 
described in Chapter 3 ha^re been incorporated iiito the 
program and any one of them can be used by changhig an index 
on the program coding card in the input data. Options are 
also available for use of both the methods of calculation 
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of streamline slope and curvature mentioned In Chapter 2, the 
option hoing exercised, again, hy change of an index in the 
coding card. Iteration scnemes are suitably modified by 
this change of index and no further chaaiges in the program are 
necessary. It is also possible to feed the streamline slopes 
and Gurvatui-'^B as input data into the program. 

A number of options are available for the input data 
required by the program dependii-ig on the nature of the data 
available. These have been described below along with the 
^at required. 

The general specification of the turbine consists of : 

1 . Number of spools. 

2. Gas constant of the u orbing fluid. 

3. Mass flow at the turbine inlet. 

4. Plo’w conditions at the inlet of the turbine (total 

temperature, tote.l pressure, and flow angle as a 

f un c t i on of r a d iu s ) . 

The spool see cif iccit ion consists of 

1 , Eotat ive speed , 

2. Po'? 7 er output if rotor exit relative angles are not 
specified, 

Finally, the spool analysis variables consist of % 

1, Number of stages. 

2. power output of each stage if rotor exit relative 
angles are not specified. 
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3 . 

4 . 

5 . 

6 . 

V. 

8 . 

9 . 

pr intec 
1 isted 

1 * 

S. 

3 . 


opecific beat at spool inlet and each blade row exit 
stat ion. 

Annulus d:eometry and axial position oi each station. 
Kass flow and total tetaperature of the coolant added 
if the turhhie is cooled. 

Streamline values of the mixing coefficient for each 
blade row if niierf ilament mixing; is specified. 

Whirl velocity or flow an-;le as a function of radius 
at es,ch stator exit. 

Streamline values of the power output distribution or 
the rotor relative exit angles. 

Stage efficiency, rotor efficiency, total-presBure.-loss 
coefficient j or hinet ic-energy - loss coefficient, 
and additional loss factors when desired, when the 
losses are not calculated from the loss correlation 
a,vailahle in the program. 

The output of the pro,^,rani consists entirely of 
data. The information included in the output can he 
a.s follows ; 

All input data. 

Tabulated streamlhie values of the flow parameters 
obtained from the converged solution. 

Tabulated streamline voAues of the mixed and/or cooled 
flow parameters for a blade row* 

Tabulated streamline values of the performance 
parameters of the stator snd rotor blade rows. 


4 . 



4:0 


5, Mass averaged perforraanoe parameters for a stage, 

6. fabulated. mass averaged perf o.rmajico parameters for 
each stage of a spool. 

?, kass averaged performance parameters for a spool. 

8» Mass averaged performance parameters for the turbine, 

9. In addition strearaline values of the flow parameters 

through each iteration of the computation can also 
be printed if desired. 

Some inbuilt checks have been incorporated into the 
program to continuously monitor the computational process, Ihe 
purpose is to abandon the calculation when it is clear that 
a converged solution may not be possible and continuation 
of the computation is not worthwhile. Such cases occur due 
to numerical instability and oscillation of the solution, 
for example, v/hen ihe design mass flov/ is close to or 
greater than the choiring mass flow at the station may lead 
to oscillations in the inner iteration loop with no convergence 
of the mass flow. Other cases may be when it is not possible 
to satisfy radial eauilibrium and the continuity equations 
simultaneously, oscillations in the outer iteration loop, 
and so on. In all cases streamline values of rhe flow 
parameters at the last iteration aro printed along with a 
message explaining the reason for abandoning the computation. 
Nine such messages are incorporated in the program and coyer 
possible problems that may occur in the computation. 
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4 «2 Inlsrf ilaLient I-Iixi 


. -A 


prCoSiace Ox Ipi^e ^rr^-dieiits of th© flow properties 
in t’QG liCT/; fi&li v^ould result in flov/ mixing that would try 
to reduce tuese gradients. To integrate the effects of mixing, 
into one comouter orogram, a simple flow model has been used 
hy Carter et, al. (l6) w hich ha-s been retained in this work, 

Ho attenot has been made to represent the oomplox processes 
by Which mixing talres place but to q.ual itatiwely represent 
the mass flow mixing ^rhich occurs betv/een the stream filaments, 

The considera.tion that if a total pressure and total 
tempera-tiire profile is modified by mixing, the new profile 
would have the original mass flow v/siahted values leads to 
the formulation of a mixing model of the form t 


Pt^ - (1 - uJ P 
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(4.1a) 


(4.1b) 


he mixing oaraneterj jijjj is spocifi-ed j-or individhal 

ti 


streamlines in the iii.ou.t '■vata., but the same voilue is used lor 
both total oressure an'? rotal temperature. 

having chosen a mathematical formulation for mixing, 
it is necessary to decide at vfnat point it is to be introduced 
into the analysis. Since mixing is related to the flow 
within a blade rov;, it would appear logical to specify mixing 
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parametero for 'blade rows rat Per than calculation stations. 

She modification from streamline to mixed values of ah solute 
total pressure and total temperature are made downstream of 
the plane at vrhich radial equilihr ium and oontinuity equations 
are sa-tisfied and the mixed values are used in the inlet 
conditions for the follou^mg hladc row in which mixing is 
assumed to .occur. 


4 . 3 Ho dlii cat ions t Cpmjout e 

Paring the initial use of the computer program, it 
was found necessary to moc’ify certain parts of the computer 
program of Career et. al, (29). These modifications are 
discussed below. 

4 .3 . 1* Im provements wne n ch_olcln^y floyi ur®. 

gnjeount ergd 

Ise of the basic cOiiiputor program of Carter, et.al . 
(29) at stations where chokiUc. occurod (usually the 

first str.tor exit) resulted in oscillations of the solution 
ana consequent ai^onaonin;, of toe- calculation. Closer 
scaruttoy sioned tuat the reason for those oscillations was 
too large a chsi«se to th® meridional velocity from one 
iteration to tiie next of the Inner loop. This resulted in 
the mass flow osolltotlny shout the oholoine mass flow without 
converyenoe. fhe obvious solution to the prohlcm was to 
use aanplnf teotoidues on the meriaional velocity 


,( 2 ) 


IQ 


yd) + (1 _ X) 


(4.2) 



43 


where value of the meridional velocity calculated 

for the new itoratioii, t'lo value for the previous 

iteration, and tho '/■alue used in the nevY iteration whicu 

is the damped value, A value of i hetwoon 0.7 and 0-8 was 
found necessary to effectively damp th<= oscillations . It 
must he emphasised that dumping is used in tuo program only 
when oscillations are sensed by it. In other circumstances 
the value is nsed 'uith no damping and thus c.uiclC5r 

convergence to tne specified mass xlo\/. 

4-3.2 ?he^ usG_ 0^ foiofi 
sp^ pif i oat ion 

The use Of the program vYhen the vjork done by the 
turbij-M is spacified rs^uires also the speciiioation of the 
distribution of the v.'orlc across the blade hoifcbt. In other 
words, the loadto,,, of eae'i section of t>- blade has to he 
sptoiflod. This was found to bj Inoonvonient booause rarely 
is the exact loadins pat torn Imoiro Cor a performanoo evalua- 
tion. However, the velocity triangles of the turbine are 

iv h fu-e rnfnr relative ex i-t anflos provide u 
usually kno-m and tnc rotor r..raorv 

^ ^ pine ororram. The equation (2.27) replaces 

better mpui^ into i-uc proo-t-t 

/r ■- +’n' T-'i n 1 v't' is and fchc stagnation 

the equation ( 2 . 26 ) rn -cim- o-uiaiy^i- a. a 

j-H -j pon-i’ nr Is calculated irem 

temperature drop across uie rOtor lo ca 

G auction (2.25). 

•ice use of either the worh done or the rotor relative 

available to the user and 

exit ancle as an niput fe an option avaiiaPlc 
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the option ie Gasil 7 oxercieod by a cliuiigo in index in the 
initie.l coding card of the inant da^.a. the /or}c epee ificat ion 
can be used for tnrbiac desi jj and the exit angle specification 
for the do 3 ign-analy s is , 




♦ ♦ 
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PJ-^ 

s^tre aisl tne^ 


PP2-PP^^ 

slppe^ arid* curr at iir:-.g. 


oyalUG-te 


'ihe use of spline fitting teclmiaues have been 
adopted by most ant'nors for calculating tjio streamline slope 
and curvature (Novak ( 9 ), frost ( 12 ), Silvester and 
Setheriiigton ( 15 ), among others). Ihe metuod has been 
incorporated in the present computer program. 


Naneshyv/ar , et . al. (30) has made a cribloal assess- 
ment of this method for use in the stroamline curvaturo 
tcehniouc. The conclusion is that fche instability ouat aiaj 
incur in the method is due to oscillations of the second 
derivc-tivo of the epline curve vdiici: is used to calculate tlie 
str'camline cur'vature, To overcome tuis problem, 

Qaneshyiiar (3o) has sug posted the une of two spline curves, 
tbs first one fitted to points giroujh '.diich the strcamlino 
passes, from which the ^--trcamliuc slopes are calculated, the 
second one fitted to the c-lculated slopes at the given 
asial stations and Is used to oaloulato tlio second deriratlTC 
and hence too sticnirellnfe ourTOture. Inc motnod, oall&d tho 
-BouDle-Spline W, is desorined in detail helon. 
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15 c radial poiiiits -throu^jh v/.i ich the streamline is to 
pass along ‘./ith their aj;:ial positions form a set of data points 
vhiicn is denoted hy A spline 

curve is made to pass through this set of data points and this 
curve csn 00 defined ajialyt ically in the interval 
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where k = 1 , 2 , 3 , . . . jh , 

21ie grad- lent at any poixit in the just fitted spline 


is comouted from 
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(4.4) 


i-jo next stfej is to replace the ordinatos 


-PI r - ^ rU^ ) r’(’- ) f’(xT,) respectively, 

hy the siopes X' [x^) , 1 ^ ^ 

A aec^ond spltie is i. en male to oass tlU'ougb taa neo set of 
points U^, f'(xp), (X2,f'(^)). •••’(^<^'(55)). S'nis 


second spline oar. bo derined in t-ne interval o. 

I-I-. + -g-- + (-T - 
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Ac:,ain, bjie {gradient at any point in thiG spline is given by : 
( - ' ' ■ ?. 


P’(x) = 

li. 


.tL 


- Q. 




(xj^ - x) ' r'(xp 


2 £ 


+ (• 




“ ^ 1^-- 1 ^ ^ 


'k "k 

k~i r 


- 0 

\ 6 '' 

(4*6) 


Ibe cnrvatLire at a typical point (x, . rj is thus 

cv iC 


cj;),, = p'(x,p/ci+ 


(4.7) 


I'he coefficients G and are evaluated by matching 
tne slopes and curvatures, of the spline, of adjoining 
intervals a,t their coumon points. This leads to the solution 
of a disy onal matrix, which is easily done on a computer. 

Two end constants have bo specifiet^ and the choice made here 
is that Used by Ka,tsanis (4), The second derivative at 
either end point is chosen one-half that at the adjacent 
point. This provides the two e:|uatlons necessarjc to complete 
the set of equations for the solution. 


4 .d 

The computer prograai ir. composed of a main routine 
and twenty kvo subroutines. A list of these subroutines along 
viith iheir main fimctions is gi'^-en below. 

1. IliF'T - to print all the input data. 

2. STTuAG - to calculate the slopt-s and curvatures of the 

hub and casing streamlines . 
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3. SPISCrlT - to determiiie the speciric heat at constant 

pressure j specific heat ratio, snd related 
parameters 8,t each station. 

4. PO,.TjR - to dereraiine the drop in ahsoluto total 

teniperrtiL.re across eacii streamline of a 
rotor v/iien tee ^orh. done hy the turhine is 
specif lei . 

5. PO hid 1 - to det ermine tne drop in ah solute total 

temperature across each streamline of a 
rotor when tne rotor relative exit an.t)lee 


6 . b iPIP 

7 . STiFAL 

8 . YTUaTJa 

9. ihiO-QL 

10. idilhi 

DUBIV 


are specified. 

to obtain the initial estimate of the radial 

position of each streamline t 

to obtain streamline values of the items 

■required for the solution of the rarlial 

e q u il ib r lum e tjuat ion , 

to obtain an initial estimate of the 

meridional velocity at the mean streamline. 

to control the lo„ic of the calculation of 

the meridional velocity distribution. 

contains i le alririthm for the Runtie-Kutta- 

G-ill method of solviii^j a firf;t-order 

or'i ino-ry d if f o rent ial enuat ton • 

to obtain a '/slue of the derivative of tne 

square of the meridional velocity with respect 

to the radial position for a specified 


11 . 
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is, ir.-iU3 

13 . -mikTE 

14 . SL'HJP 

15. OU'i’PTJ’I! 

16 . PLO 

17. 

18 . LCITV 

19. SPLIIsS 

20. I1-4P1 


Yelociti? aiid radial position, 

to obtain a new estruate of the meridional 
velocity a,t the mean streamline, 
to obtain streanQ-ine values for the quantities 
tabulated in the output which have not already 
been ob t aine d , 

to obtain streavalitie values of quantities 
for calculations at the next station. It 
also calculates the mass averaged values which 
a,re to be printed in the output, 
to print the resulto of the calculation, 
evaluates the total-precBUre-loss-coef ficient 
for each streamline using the loss model of 
Carter j et , al. (l6), 

evaluates the to tab-pressure-loos coefficient 
for each strcmlhie using the loss model 
developed in this vorh. 

estimate! the total-pressure-lOBS coefficient 
from the kinet tc-energy-loss coefficient 
:,'hen they -re specified, 

evaluate! streamline slopes and curvatures 
by fitting spline curves through the 
estimated streamline radial positions, 
is to perform parabolic interpolation of a 
tabulated function of one variable, if 
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21. SLOPIJ 


22. si-j:q 


p.^rabolic interpolation is not possible, 
linear interpolation or extrapolation of 
a siiii^le value is peri'orraed , 

— to oPtain tlie derivative of a tabulated 
fimction v/ith respect to an independent 
variable ea-cb tabular entry. 

_ sola!. ion of a set of siaiultaneous linear 
algebraic oquat iono by the G-ausc method, 
vritli pivot llig. 

An over-^all ilov' diap;rara for the prograir is given in 
figure 4 . It can be seen that tbe over-all control of the 
calculation procedure is maintained by tbe main routine 
;vnile subroutine HA-OSai maintains control over tbe calculation 
procedure for the meridional velocity distribution. 



applic^ko:-s jioa'ssiCi 


The coTouter jrogrcum wa,6 used to the predict 
of severs,! single axid maltistace tarhinea , Ibe 
these applications are presented in this chapter. 


the performance 
results of 


to dehag tne ooTg.ater program the results presented hy 
Carter, et » al.(29) were duplicated. This cleared, any 
poBsihle errors in the basic comnuter program. Necessary 
'modifications wc-re then made to the program in steps and the 
modifications debugged at each stage. 


■I‘o study the effects of the radial cilsturbution of losses 
on the solution, efficiencies vvere prescribed in the initial 
program runs. rliis provided Oi useful insight into the effects 
of the loss distribution on the solution, Ihe first of the 
problems encouait&red was in the prescription of the work 
distribution aloir; the blade heigut. The turbines under 
consideration were designed with .. constant loss all along the 
blade height. OonseQuently when a radial distribution of the 
loss vms prescribed, the design worlc distribution wea no longer 
valid. fhie n-as immediately evident in initial program runs 
by an inability to find a valid solution to the radial equi-- 
librium equatio'-i in the hub and tip regions due to overloading 
of these sections. Since the efficiencies prescribed in the 
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Ijub'aud tip regions nere loser tlion the average efficiency 

in the design of the turbine, to obtain the design work 
laiger turning of the flow at these sections was required J 


larger turning of the flow at these sections was required and 
consequent overloaairn . lo oo^ercorae this problem required 
a';redistribution of the v/ork lone along the blade height till 
the rotor exit flow angles closely matched the design angles, 
ms proved a cumbersome procedure. The use of the rotor 
ixit flmv angle instead of the work distribution as an input 
to the program was an obvious solution and a subroutine was 

added to facilitate tnis option. 

me use of radial loss distributions introduced large 
^^dients in the flow properties which increased as the 

^calculation proceeded dov/nstream till their profiles 

^ n -aa nn longer possible. However, 

‘"so distarted that a solution -.vas no longer i 

I' on-u- in real flows. Mixing of 

' such a situation would nox oc.u. 

„ /^n+ ihi crradient . The 
W"' 1 j. T 4 - i^’T’cen oux bni-s o-^ 

i the flew would teiis to ..^raro*- 

,-.c.pribea in Chapter 4 vms used to 

; simple Eiixin& mcu. — 

• It WPS found that pres- 

..-simulate the effects of mizcuf. 

M-,!-..- in the floc' for each otage produced 
cribmg 50/i- ^ • oi i 

„ of nixing was maintained in all 

best results. i-aini i- 

. .a,, ., 3 sults OX wbio'n are presented here. 

" program runs, ua-*- 

■ 4 .1 2es^'t__pa§e_l 

. -13 -tame turhine with large huh-to^tip fadi^ 

A single oOc5.ge u x rST/A 

Urst test case (Ref ^ 


rs-tio %?as used a-s ^^le i 


A .R31l!S3 
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details of this turbine are presented in Figure 5, fhe 
progra,m was run \^ith a prescribed work distribution, 
prescribea rocor rel^-tiye exit andies for tho present loss 
model and the loco modal of Ref, (i6) . 

Figures (6) arw'! (7) pres nt the r-adial distribution of 
outlet flo7i angle and stagnation pressure distribution. Both 
these distributions are predicted v/ell by both the loss 

models. I'his indicates that the effect of the loss distri- 

« 

button is less important in the case of large hub-to-tip 
radius ratio turbines, i*e., turbines with small blade heights. 
Figure (8) compares the efficiency distribution of the two 
loss models. The figure clearly shovvs the di If erence be'’cween 
their loss distributions. The present loss model shows a 
maxiraux'ii efficiency around the mean line as would, be expected 
in an axial turbine \/hile the loss raodel of Ref . ( i6) shov/s 
a maxinura efficiency at the blade tip which is unlikely to 
occur in such turbines. 


Table 1 couuaros the prodictod o-vorall performance 
parameteirs with expo rime nt ally measured values. The loss 
■model of Ref, (i6) predicts e. higher efCiciency and lovrer 
pressure ratio vihile the present loss model predicts a lower 
efficiency and more accurate pressure ratio. Better 
prediction of the \7a.ll static pressure is also indicated by 
the present loss model. 



cperfe-tal Present loss Present loss Loss model 
model (’vvork model (angle of Ref. ( 16) 
spccifica— specif ica— {v^fork spe— 
tion) ci float ion) 


■ total 

p: 4 Xf ii? i C3.1 cy 

0.92 

0.89 

0.89 

0.93 

. •) .'lEiptaJ-to-So'-al 
-pressure ratio 

1.909 

1.897 

1.909 

1.84 7 

stat ic 

0 . d 'i5 

0 .462 

0 .463 

0.481 

■' static 

'-/v- pressure 

C .44: 5 

0.461 

0.462 

0.480 


__ , _4 ,2 l!e 0 t Ca se 2 

) A kvo stage turLine designed for a low cost turbofan 

used for tbis test case (Ref. 32). Design details of 
t-be turbine are presented in Rigure (9). Lhe program was run 

i * ■■ 

-'Pptising a specified work distribution with the present loss model 
S-band specified relative exit angles for the loss model of Ref.(l6) 

, 2^-.' fhe redial distribution of the outlet flow angles at first 

atnd decoiid stage exits predicted by tne two loss models are com- 
- V pared vith experineBtally measarad ralues (Sigures (lO) and (ll) 
respecti'rely) . It is evident tnat the present loss model shows an 
Improved prediction of the trend of the first stage outlet flow 
’1 angle, large deviations of the predicted outlet flm angle at 
i’-' the second stage exit, in the hub and tip regions, from the 
)';■ experimentally measured values are seen In Pigure (ll). However, 
aef. (.52) points out that the measurements showed a large amount 
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of underturnin,3 in ihese regions. A large amount of turning 
at tnes6 sections would have resultod in o. distribution close 
to that predicted by the preoeiU loss model. 

Table 2 compares 0 Tera.ll performance parameters predict- 
ed by the b-'o loss models -ith ernerimentolly measured values. 

The loss model of bef. (ic) jredicto a very Aovi total pressure 
rat io . 

Tcbic' 2 


Szperimental Loss of model Present loss 

of Ref , (i6) model 

Turbine total 


efficiency 

0.93 

0 .9073 

0.914 

Turbine static 
eff ioierxcy 

0.00 


0 . 777 

Tot al-t o-t ot al 
pressure ratio 

3.75 

2.525 

3.51 

Tot al-t 0 -s tat i c 
pres.sure ratio 

4.6 

3.845 

4 * 63 

Stage one total 
efficiency 

0.93 

0.9 5 

0.908 

Stage two total 
efficiency 

0.91 

0.33 

0.903 


2 T est C ays e _ ^ 

A three stage turtinc vvith very high stage loading 
was used for the third test case. Resign details of this 
turbine are presented in Ref. (33) and tue results of 
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experimental investigation in Eef , (34) . Some of the design 
details are presented in 5'Lgnre(lE). ihe program was run v/ith 
specified rotor ral?.ttve exit aji las. 


ihe predicted radial distr ihut ions of the 


XU bciJJ— UU- 


V/ 


Total efficiency for the firsc, second^ and third st£ige exits 
are compared v. ith experimental results in Figures (l3), (l4) 
and (15) respectively. Good agreement is shown in all the 
cases. The de /-iation of the efficiency level at the mean 
line suggests that the estiaiated profile loss using the Ainley- 
Mathieson correlation is too high (Horlock (2o)). The predicted 
'/all static pressure distributions compare very well with the 
experimental results (Figure (l6)). The predicted outlet 
flow angle distribution (Fi.uire (17)) shows deviations from 
the experiment ad- ly measured values in the hub and tip regions. 
Deviations of the same order are also seen in the results 
preseiited h/ ilirsch (V) and larger do'^/iations in the results 
of Renaudin and 3omm (l4), fcr other turbines. 

Attempts a,t obtain'.: a joint ion using the loss model 
of Ref. (lb) faileu in * nis tesij case • ihe apparent reason 
v/as that with niglily loaded, turbines ^ the solution is 
sensitive to t„’s radial dis iirioution of 

Tabic 3 compares tne predicted overall performance 
parameters vritht oxperhflcnta,! results, Good agreement is seen 
v,'itn the experimental results. 
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,ble 



Experimental value 

Predicted 

Turbine total 
efficiency 

0 ,8ci6 

0.867 

Turbine word output 

3Cir,03 lip 

29 52.81 dP 

Total pressure ratio 

3 . 43 

3.42 

Total-to-gtat ic 
pressure ratio 

3,843 

3.839 


Test Oaa e 4 

A four stage turbine vith very high stage loadiiig was 
used for the fourth test case. Design details of this 
turbine are presented in Ref. (35) and the results of experi- 
mental invest igcit ions on It in Kef. (36), Some of the design 
details are presented hero in figure (l8). The progroca was 
run wit’.i the rotor relative exit angles specified. 

Figure (i9) compares the radial distribution of the 
Total-to-Total efiiciency v/ttli the experimental measuromonts , 
While the trend o the experimentally measured distrihuuion 
has been predicted v/ell, the overall eCficiency level predicted 
le lov/er by about , Onco oS-onni, tin difference in 
0 y ( 23 ^ 0(3 3 /b t ii 0 mean line 'uould indicate a too uigh 
profile loss predicted by the prcsint loss model. Good 
prediction in the trends of tlie radial distribution of the 
stagnation pressure and outlet flow angle are seen from 
Figures (20) end (2l) respectively. 


file lov/er stagnation 
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pre,3sare could "be duo to the iowei: predicted efficiencies. 

Tae deviations in tijc outlet .flo ' an, pie are of the eajiie order 
as i.he results snoem t-y nivseh (7) for other t nrh ino s . Good 
predict Lon of tcc 'vail static prosouro distribution ic seen 


from Figure (fS). 

fable 4 cofipares cho orodicted overall performen.ee 
parameters -:ith e^cperiruuntal results* As mentioned earlier, 
the predicted overall efficiency is lov/or than the meacurGd 
value , * 



Table 1 

Exp 0 r ime ntal 

Predicted 

Turbine total 

efficiency 

0.853 

0.835 

Turbine work outp 

ut 3400.00 HP 

3380.10 HP 

T ot al-t o-To ta.1 

1 


pressure ratro 

2*66 

2.875 

T 0 1 al_ 1 0 - s t iit i c 
pro sou re reitio 

2.9 55 

3.221 


Results obts,ine-| from irf tiar program runs using a 
preacribod offioiGricy distribution arc prosented in Fipnros 
(23) and (Sd). Good agrocfn.nt ’vith experimental results 
in all cases rjhoved the iiiiportanco of the loss distribution 
on the predictions. In £i hid to study the effect of a change 
in the prescribed mixing pattern, Figure (2 5) sbov/s the effect 
on the static pressure diotrihat ion of prescribing the mixing 



58 


coefficients by t'lo jnetliofs. Very little iiffcrcnco is 
observed by prescribing mixing l'g etage exits and proscri- 
bing 3cy^ raixiae at all blaie rO'j exists. 

Attorapts at obtaining solutions using the loss model 
of Ref. (16) failod in t lis tost case also. 


lo substantiate tbo stat..-ment that the doviation in 
tne rosulto is due to tiio high profile loss from the Ainloy- 
liathioson correla/iionj the pro; ruixi ’•/'as run after reducing the 
profile losses commuted ■'uith the acov'-^ correlation by dOfo . 
Secondary losses weie kept uiichanged figures (^i6) and (27) 
sho;/ the nevr radial distributions of total pressure erul outlet 
flow angle coanared with experimental measurements, Oonsidcra- 
bl 3 improvement- in the predictions 'are seen as a result of the 
above mo 1 if i cat ion in tea loss correlation* 



C-liAP'ljjlR 6 


OpRCIiyb^IfS 


An attefiot UaG been xn tbe present invest at loiiS to 
evaluate CLlflerent lose models for axial turbines in a 
streamline curvature comijutino ssbene using ttie lull 
radi&.l en/uilibrium nation, lour single and raultlstag 
turbines ba.ve been analjsed. lue design details and 
measured data 're re available for these turbines. 

Tbe loss model proposed in tb is vorlc to take into 
account tbe radial distribution of losses in the airaulu 
has been very successful. Ihe predicted v/all static 
pressure distributions and the radial distr ihutions of 
stagnation pressure ajid efficiency agree very well with 


the raeasurecl ones. 


The loss model of Garter jOt. al (l6) failed to produee 
solutions ior highly loadjd lauit istagv-: turbines. This 
may be attril'uted to 'i r-i ael distribution of losses 
quite d if l'.^ r.'.nt from uuosc oh'tciiiio;.. in praotiGC, fj-iis 
model also prodictv an o/erall efficiency riiich is 
higher tbiui tao moor.uaed one. 

It v>' 3 r observed, that tae Atnley-iViathiesoii correlation 
predicted too hig’n a value of profile loss for modern 


turbiaics. Infs load to a prediction of a lower overall 
efficiency for the test 'burhiiios. A reduction in tae 


level of the profile loss p.roduced conslderablo 
iiflpr ovemcn'bs m tao results . 
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The eva,luc/t ion of streamline slopes and curvatures 

from the c-xlcL cr.s in^; slopes end curvatures was found 

satisfactory In vljn of tlxu food pr'Cdictions obtained 
for mult is t ar , e t urb ine s * 

Several modii ications liave been made to tde computer 
proijraiii of Ga.rt^ryet. a-1 (29) • I'O'-' addition of rotor 

relative crit ajn_,le subx’o Lit trie lias made it cuitablo for 
both anal;yois and design computations. 





;^L A'- .^Pl ^ Ap. ®iiPiAAAA^P*P^ 


Tlie cont ii3ui'*'r erLia+ 'oii c"-i 'je v/ritterj as : 


li P) 


+ V . i 


P Ds 

In cylindjL’ico.1 oo-o id mates. 


1 8(r-Jj.) 8;/^ e;.;_ 

V . ‘-"it,'''’' + "''t't" *r . 

r dr rSO fix 


X 


(Al) 


(A2) 


Also, the rlsrivative of a qu&Jiti by in tlie direct ion of flow 
(tide s-direction) at atiy mstant in time may be expressed by 


§iJ. 

5s 


9lJ_ + 1^. ^LX + !ls . lLI 

dr V/ rfie ’W fix 


W 




(AS) 


or ro-arraii;^:ui,:';, 

iL 1 A 


A.i _ .'j; . ^-Ll _ 1 . ?Xi 

' fis fir '*' 


6 


W ' rfie 

X 


(A4) 


036 equation (A4) in (A2) to eliiulnsito t;.e x dorivatlYe 

■jVe -^d^et 


1{ 


r * b: 


rib 


(p-) 


fi 

rfiO 


v-,r / 

I i 

X 




+ r 


X 


’’ fix 

(A5) 


Since 


vr 


tan <j!i aj' 


v/ 

01 

i.C, ' 
X 


tan p 


e;j nation (nS) can be rev/ritteii as 
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^ . 7 / = 


-L't. lo 'bha^t 


PS,^L ^ J j. L 't an g 

as- iae ■ 


• dx 


(16) 


2.U. = MJ. + -„, Si Jl 

^ O \J OS 

l 1 gIi 6 iciriG'CQr-„‘t ic x'GXti't; ion 
can )je v/ritten 

ai-i = 2L 1 1 ai J. 

53: Ds ¥ * ot 


(17) 

, e<|uatIoli (AV) 


(AB) 


^(^■uauion (A8) caii oe '-.sed '"or tho lasx iGrui of Gq_ nation (A6) . 
Also since 


= W cos 


(A 9 ) 


V'/e .'ret , 


Ac D 

r = TS“ (<: cos 


7 • af 


(AlO) 


17 e a.lso laiov/ that 


W -U- ^ iY 

Its al Di'l 


(All) 


Ihus, 


do ~ ;/ 


D cOo 5 
jlLl 


7 T 

D'j . 

Ill 

ii 

A"' • 

Dm 

(A6) 5 

we 



1 z 
J 6t 


(Ai2) 


Suhst itnte eooation (Al 2 ) in (A6)5 we i,et 
y 1 ^ w + 


\Y et 

X 


(A 13 ) 
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J’'or an steady Xlov; 

Y ^ y ^ X Jl). 


Dn 

ec_^ ■ 


+ '.V 4- 

r r JJiii 

lii 


(A14) 


X 2? dj? 

We next aopj-y t]ie '101061141112 equation iu tjie Clow direction 

6p 3 2 

- r * = TvW (. 0 “; - 1 


p 6 0 lit 


a) rw 


(A15) 


equci't : 

'■ on s 

(AS) 

nod 

(^11) 

Dp + 

1 

_ 

5£ _ 

r*> 

1!“^ 


De 

p 

dt 

' m 

Dili 


DW 

I 'V ^ ■ 
u Dx 


i"- -- rW 


(A16) 


[Pbis oon be rewritten as 


,1 

p 


- D;7 y 

IHE w 

* Ss' “ ' m Dm"^ “ r r 


p 


u dp 1 &£ 

’ rde P &t 


(Ai?) 


Tor an axis ysnie trie, steady ilow 


■p ' Ds “ IU 


JJil 


w 


V 


— i7 ““ 

Din r r 


Also. 


Dpps 


(Ai8) 


(Al9) 


ibis relation alno - "itb cpnatioxi (Al) , (Ai 4 ) , and. (Ai8) can 
be used tc obto.m an express ior for Dl^iia tliat does not 
contain V, Dp/Ds or Op/Ds. ie pet 


1 ”-r 


V 


S 

-r — 

" TiTrl 


d ;7 


"X 6(Eta£,j} 


DA 

se£ £ , ^ _IS. 


(ASO) 
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Using the definition of Maob nmi'ber, ^'^0 “* 

in er^uation (A20) and rearranging, we get 


(1 ~ u^) ^ - w 


m 


‘M 




sin 


, t. dii 4> , a in 4 ^ 1 6 4> 

^ ^ v- oos dr 


m 


(A21) 


Since V , equation (ASl) can be rewritten as 

I'li iU. 

o 


t i 


m 


7. 


m 


;.l 


(1 + H? + .v^- + taa-^ ■}> 

6 eop (() ■' r or 


Dm 


-'ti’ 


(1 - Mp 


(A22) 


sin 4 > 



■Apj2riid ix: 'B 






when the total jjover output of o sta^^e Is specified, 
the total tci^iperatare distribution at a stage exit \/ill have 
to satisfy both the apeoiticd j0\/3i-' output o,acl its distribution 
across the annulus, oince aiiitiallj'' the distribution of 
inass flov/ throughout the aiuiulus is unlcno'.ni until toe distri- 
bution of f'Jeridlonal velocity has been erlablished, the power 
distribution is specified by nondiiionsionul power fiLnctions 
versus tne nond imens ioiial mass flow function, v/(r), defined 

271 f V,. r dr 

r ,A. 

(^23) 


w(r) 


w, 


If tliG tots.l pov/er out out Bpecified is liP^ (Uorsepovfor) , 
the total temperature Iron a I ^ throiign the rotor must saticfi 
the eoyiation. 


r f 

IIP r. 5,692 5 c / .bs 

^ ^ o 


(A24) 


Kortnalisln. . ecamtion (a.M) with rorjoect to the i;otal po’uer 
and the total mass flaw, lee-ls to a defjn.it ion of the power 
function P(Yr(r)) expres.:od as 

v/(r ) 

p(Yr(r)) = 5.692 5 f A d wir) 


0 


(A25) 



Dixi'ertint ul‘c mg, equauion (-/E5) 
dimensiona.! mss function 

total te^iipe nature drop, j’or tlic 

■tiPrp r 

0^3 0,09-0 v.’m c, L 

j * 


■;itli respect to tde non- 
yiclcir- an expression for tlie 
2 1 li r ; t r e aul in e , 




(.•i 26 ) 


Tile power function mrsas 'fno mass flow rmetion v;ill le a 
"basic specification for pov/er distrilutioii Crom which the 
total tenipercvture drops are o"btalned, 

\/hen tile rotor relative exit angles are specified, the 
total temperature drops are o'btained from the Euler work 
equat ion 



Appgi -d ir Q 

CjDe^f oi\ t rig _ J- if -Cp ir.nt la l e i. | uai iona 

io siriolif;' the lo_ic oi? jL.-j coaiguker nrosraii, a 
standard oroceduro has hesn adopted tor the colution of the 
flov/ fiold at each inter hlade ron station. The different 
types of stations and tag varior.s optional specif i oat ions ane 
tahen ii^to account hy iuodificat ions to the tvelvc ooefficientn 
appeanin{£, in the three difiorcntial e\.|U.ationB (ilquat Lons 
2.30, 2.31, and 2.32), This appendix presents those 

coefficients , 

P ir s t st at ion inle t 




- 2c 'T ^ 
uo no c 0 0 


COL. 

.*1 0 


+ (y- + v:j p 

r UU luO a. 


■ IilO 

+ y Sllr <})„ 

DIO 0 Dri 


ineciried function of radius 


The total preccarc -ouLd oe a specii 

hence, the coefficient, of enuation ( 2 . 3 l) one 


°gi = 0 


1.0 



68 


trn p coc 


Ogg - 0 


Cgg ~ 1,0 


= V 

34 mo 


Statoi' exit 


cos 4 > da 

2 ^ di’ 
cos 6 

0 


d 0 

tan a sill A 

•^0 ^0 dr 


=: 1.0 


0i2 = [\u + y,!i - «o 


1' 

P 01 


2V^, oos 2Y; 


' A -4 **> V j I j * J f't I I J. 

^ ^ .r2 N 1 ;j-oi 


+ {¥ + V'" ) 

^ U1 ^ ‘ml-^ 


+ j s ij'i 


1 Dm 


Assuming, t.dct t 'le deriv'^tlve -™-“ cax ijc-' oxcresscd as 

dr 

dY.: ::V 

dr " "II dr ^ 4'3 dr ^Y4 


tile coefiioionts oi' equation (2.3l) Dccoujc 

^ ^1_ 


(Aj-) 

■-0 0 __ 0 1_ _ 
^^“2¥ T „T 


Pm Pi V 

+ (:i'-’^-‘)(l - 

Pno Pm 
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C 


22 


1.0 


^23 ~ 


V (foi) ( 
Ul ^ ^ 


oo 

^ ~Y s' 


V 

r ^ 


■ol 


0 


01 


q’"'") 6!^g 

■" o 1 


0 . 



do 
- 00 

24- 

P 

„ 

OO 


°31 


ta]i 

■b;9« 4. R, 

^ ^ OOE 




ml 

n 

"32 

= 0 



"33 

= 1. 

-0 




r 

cos 

S4 ^ 

= tai 





cos § 


J. Pp _0 1 

:.'b^ rfr ■" •‘' 

jJ* - 

ol 


/ 2 p 4, <11 , 

(4 , + ) -4^ ^ri 

Lil nil'' S , dr 


01 


- ( 1 - ll- 


'P 


00 


Ol 




dr 


- tan sijj 




i!j 

dr 


Stage e:x it 


°11 = 1-0 


'^12 “ (yTI 


tl2 + 4 


■^'^p -'02 


C 


1 




G 


14 


%2 


111 


Qy^ 

V o p n 

-F- ^ + 4) f 


o2 


(IS 


DV, 


+ Y. 


mg 


in (j) 


mg 


2 ]>ii 



?0 


dY 

Again, assuming that tne derivative 


oan he expressed a,' 


/i 

dx 


R 


d ? 


dr 


- 


2 

ni 2 


dV 


u 2 


+ 0 . 


T-1 dr " .3 .lr' ^Yd 

tlie coeffi clients of oql.., li-'n (2, 31) are 


G 


P ’-. 


__0 2 ^ 
21 " ^ 2 E ' T , 


^ P^P Po 

- + ( A ^) (1 - - f -) a 

^ oSb ^ o 2 ^ 


n 

-j 


c 


23 


1.0 

P'o 

"■" 32 s _1 

B d 


U . 


(„ 0 £ L ) y , 1 

^ o 2 o 


■R 


■R ' T ' 
-^02 




+ (;^-) (1 - “' r “) C 


^ o 2 s 


Po 2 ' 


' Y 3 


0 


1 




1 


24 


Pil dr - 2R T-g 


2 Uc 3 ^ w 

^ X 


I ^ 2 \ 1 

(u^ - u ) 

1 ^ 0 l 


^ ol 


Pi 2 

+ — f - *- 

VJ A -* )J 

^ 2 E X '); 

o 2 


(u 


2 

2 
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" o2 


Yp 

io 2 s 


2 R T 


- (^) (1 

Po 2 s 


R) , 



Lip 

1 

Y, .,) 

Ml 

lX<f3 

■“ og 

9 

2 

OCn 

nlic. 

\2 

Po , 

" “r ) Cy . 

Pig 


^'^ u 2 ^ 


a ' 

"' o 2 

dF 




d'i 


o; 


_ o2 
dr 





',!ixen "bje .orK oatput of tlie turijine is spccifieri^ is 

iii effec-'i a 1 soov;k function of rai Lus . ThiiSj 


n „ 

•-' J. 



^ c 


G^-- - 1.0 



V.1921 t!'e rotor reloTivo e^^it an^’les are spocifred, the 
co-ifxi-cisnt s ox er^iiatioii (;;.32) are similcir to those at the 
star or e^d-t. 



tan 


2 


c,,p = 0 


%'i = 


Si = O2 


coi 


I 

'2 


^2 


dr 


d 4 > 

tan ti* atn tg 


+ LI 


'=■ c,os 



oend i3’_D 

S®, cj'^e •: olut ioo of ord mary 

dj-f/erentialj^lkfci^ ‘ 


I'D 3 .1 iffer-atiai e.'iuatioii io ox tlu; 


1 orm 


dT/ ;r 

ul' 


:(i, TJ 


v/here fCj-' i's obta iii.a-1 xvc;,i tlio s iiQaltaii£'OU& solution of 
t u e r a o ia 1 e q u iJ- ib r i urn e q uiil i on and / o f ) ub s id ia r y d if f fe r en t i u 1 
eqUat loi'iB , 

Gi;,'-en ino vR-lue of ibo meridional volocity nt one 

jt re oialirie , ('/.)., the unknoY/n value at tiie adjaoont stream- 
mi 2 

linej (V •)'-j is deteriained in four* stages ; 

TU-L j.^ 
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(3) 


^ ^ -fc) (1^3 - ’ip) 
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i-ij? {r + (y > ^ 

X.i ' ? ' yi*i/’i_OJ 

^ --tj- *i^ro 

'1 + 3 I ( 1 - ^v) ty - q^) 


- f’'ai\-3 + S' - -(Ig) 


" ^-3 - 
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